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■ ABSTRACT 

^ : 

^ |. Context. There is growing evidence for the presence of blueshifted Fe K absorption lines in a number of radio-quiet 

I ' AGNs and QSOs. These may be fundamental to probe flow dynamics near supermassive black holes. 

O ' Aims. Here we aim at verifying and better characterising the existence of such Fe K absorption at ~8-10 keV in the 

, luminous Seyfert 1 galaxy Mrk 509, one of the most promising target for these studies. 

C/3 , Methods. We present a comprehensive spectral analysis of the six XMM-Newton observations of the source (for a total 

, ^ , ' of ~200 ks) , focusing on the detailed and systematic search for absorption features in the high-energy data. 

' Results. We detect several absorption features at rest-frame energies ~8-8.5 keV and ~9.7 keV. The lines are consistent 

I , with being produced by H-like iron Ka and K/3 shell absorptions associated with an outflow with mildly relativistic 

^ ■ velocity of ~0. 14-0.2 c. The lines are found to be variable in energy and, marginally in intensity, implying that variations 

QQ ' in either the column density, geometry and/or ionization structure of the outflow are common in this source. 

Key words. Galaxies: active - X-rays: galaxies - Galaxies: individual: Mrk 509 

■ 

(N ; 

^ 1. Introduction et al. 2009, Giustini ct al., in prep.). Variability studies 

1 ^, . „. . , r , ri 1 1 -r 1 could, lu principle, allow us to estimate the location and 

^ . There is compelhng evidence for the presence of blueshifted .^^^ the absorber (Nicastro et al. 1999, Krongold 

O . Fe absorption hues at rest-frame energies of 7-10 keV m 2007, Risaliti et al. 2005, Reeves et al. 2008), thereby 

> ' *^o^ 3 'P"" '""^^'^ radio-quiet (RQ) AGNs and ^^^^^ stringent limits to the remaining large (order 

- QSOs. These are c^only interpreted as due to resonant ^.agnitude) uncertainties in the total mass, energy and 

^ absorption from FeXXV and/or FeXXVI ions associated to j^on^entum that is driven by these outflows (Elvis 2006). 
_ a zone of circumnuclear gas with log 4 ~ 2-4 erg s cm and 

column density ~ ig^^"^** cm~^ (see Cappi 2006 and refer- Overall, these studies are fundamental to probe and 

ences therein), i.e. significantly higher than the absorption characterise the most extreme (in terms of velocity, ioniza- 

layers commonly associated to warm absorbers (Crenshaw tion state, mass outflow rate) phases of outflows, so helping 

et al. 2003, Blustin et al. 2005). in quantifying the potential feedback impact of AGNs onto 

The velocities of these absorbers being often quite large, their host galaxies and the IGM. 

reaching sometimes mildly relativistic speeds (up to 0.2-0.4 / nr^„.\ ■ r 1 1 • 1 

c), suggest we might be probing the extreme component of a , , f .V-if '''' 2 1?^ brightest (2- 

wind/outflow that may well carry large amount of mass into flux of -^2~5x 10 "erg cm ^ s i) and most lumi- 

the interstellar medium (ISM) and/or intergalactic medium J^^!^^ (L2-iokcV -^1-3x10 erg/s) Seyfert 1 galaxy known. 

(IGM) at a rate comparable to, or even larger than the interesting science case for these studies because a 

Eddington accretion rate (King and Pounds 2003, Reeves Previous study by Dadma et al. (2005) hissed on BeppoS AX 

et al. 2003, Chartas et al. 2003, Pounds & Page 2006). niuhiple observations and on two ear y XMM-Are^i;ton ob- 

Most interestingly, variability of the absorption on servations, found a transient and blueshifted absorption line 

timescales down to a few ks has been found among some ^* f rest-frame energy of ^ 8.2 keV. In addition, a transient, 

sources (Risaliti et al. 2005, Dasgupta et al. 2005, Braito fedshifted, absorption hne at a rest-frame energy of ~5.4 

et al. 2007, Turner et al. 2007, Reeves et al. 2008, Saez keV was found. Motivated by these resuhs, additional mon- 

itormg of the source with XMM-Newton was then proposed 

Send offprint requests to: M. Cappi and eventually performed. Here we present the results on 

e-mail: cappiOiasf bo . inaf . it the systematic search and characterization of Fe absorp- 
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Table 1. Mrk 509 XMM-Newton Observations Log 



N" 


Obs. ID 


Date 


Expo 
(ks) 


Counts 
(2-10 kcV) 
(xlO-*) 


1 


0130720101 


2000 Oct 25-25 


20.7 


7 


2 


0130720201 


2001 Apr 20-20 


23.1 


10 


3 


0306090101 


2005 Oct 16-16 


2.7 


1.1 


4 


0306090201 


2005 Oct 18-19 


59.8 


25 


5 


0306090301 


2005 Oct 20-20 


32.4 


14 


6 


0306090401 


2006 Apr 25-26 


48.6 


23 



Number, Observation ID, Observation date, net EPIC pn exposure time, 
and total 2-10 keV counts for the six XMM-Newton observations of 
Mrk 509. 



tion lines at high energies obtained from a comprehensive 
analysis on the full set of six XMM-Newton observations^. 

2. Data Analysis and Results 

Mrk 509 was observed by the XMM-Newton EPIC instru- 
ments on six separate occasions, for a total good exposure 
time of ~280 ks (sec Tabic 1). The observations were per- 
formed once in October 2000, once in April 2001, three 
times in October 2005, and once in April 2006. The only 
contiguous observations were those performed in 2005. The 
data were reduced using the standard software SAS v. 8 (de 
la Calle & Loiseau N., 2008) and the analysis was carried 
out using the HEASoft v. 6.5 package^. After filtering for 
times of high background rate, the useful exposure times, 
corrected for the live time fraction (up to 0.7 for the pn 
in small window mode), were between 2.7 and 60 ks per 
observation (see Table 1). The EPIC pn and MOS cam- 
era were all operated in the "Small Window" mode with 
the thin filter applied for all observations. We checked that 
pile-up was not significant in the pn data at this level of 
source flux, even when considering both single and double 
events. The pile-up fraction for the MOS data was ^15% 
for all observations. Also in this case, this will not affect 
the results on the narrow absorption lines presented in this 
work. Source counts were extracted from a circular region 
of radius 40 arcsec, while the background counts were ex- 
tracted from a nearby source-free area of the detector of 
the same size. Given the significantly lower effective area of 
the MOS detector between 6-10 keV, and for sake of clar- 
ity and simplicity, we only report here the results obtained 
from the pn data and show the MOS data for consistency 
checks only. We decide also not to report here on the ob- 
servation n. 3 because its very short exposure (^2.7 ks; see 
Table 1) hampered any significant and interesting result. 

2.1. Light-curves 

Analysis of the soft and hard X-ray light curves from the six 
observations indicates only moderate (~20-30%) flux and 
spectral variability on both short (< few tens ks) and long 
(^ycars) time-scales (see lightcurve in Figure 1 of Ponti et 
al. 2009). A more detailed time-resolved spectral analysis 
focused on the low-energy warm absorber as measured by 

^ Results on the Fe emission line properties, for the same 
dataset, were presented in Ponti ot al. 2009. 
^ http://heasarc.gsfc.nasa.gov/docs/software/lheasoft/ 



the RGS is deferred to a specific paper (Detmers et al., in 
preparation). 

2.2. X-ray Continuum and Fe Emission Lines 

Given the weakly varying continuum, we extracted the 
mean pn spectra for each of the six observations, grouping 
the data to a minimum of 25 counts per channel to apply 
the minimization statistics. Only the data in the 3.5- 
10.5 keV interval were considered, and the signal-to-noise 
ratios, even at 10 keV, were always greater than 12 (pn) and 
6 (MOS). The spectra were first fitted with a single power- 
law model plus a column density with value fixed at the 
Galactic one (A^Hgai — 4 x 10^° cm~^. Dickey & Lockman 
1990). We checked that the addition of a further neutral ab- 
sorption component intrinsic to the source was not required 
by the data. Wc obtained good fits of the continuum with 
a typical power-law photon index of F ~ 1.55 — 1.72, but 
residuals of the data (Fig. 1, top panel) clearly showed the 
presence of a neutral Fe Ka emission line at the rest frame 
energy of ~ 6.4 keV. We therefore added a Gaussian emis- 
sion line to the model. The line was slightly resolved (see 
also Ponti et al. 2009), and we fixed the width a to its 
best-fit value of ~100 eV. Best-fit values with this simple 
model are reported in Table 2, with errors quoted at the 
90% confidence level. The fit improvement was of Ax^ > 42 
and the equivalent width of the emission line of the order 
of ~60-80 eV for all the observations. 

Table 2. 3.5-10.5 keV continuum and Fe emission line param- 
eters 

N° T Energy^ EW^ Flux^ x^/f Ax"" 

(keV) (eV) 

1 1.55 ± 0.03 6.37 ±0.03 82 ± 40 O 924/850 W 
1 1.60 ±0.03 6.41 ±0.04 60 ± 30 3.8 992/991 42 

4 1.70 ± 0.02 6.41 ± 0.02 65 ± 20 3.5 1151/1194 130 

5 1.72 ±0.03 6.41 ±0.03 62 ± 28 3.6 990/1012 67 

6 1.69 ±0.02 6.43 ±0.03 58 ± 20 4.2 1175/1189 101 



(1) Emission line rcst-framc energy centroid, in units of keV; (2) Emission 
line rest-frame equivalent width, in units of eV. The width of the line was 
fixed to (T— 0.1 keV (see text for details); (3) 2—10 keV flux in units of 
10"" erg s"^ cm"^ 



In a few observations (namely observations 4, 5 and 6) 
we found additional weak, but significant, residuals blue- 
ward of the 6.4 keV due to presumably Fe K/3 emission 
plus higher ionization iron emission components. Indeed, 
using this same dataset. together with a ^^100 ks Suzaku 
observation that was partially (for ~25 ks) overlapping the 
observation n. 6, Ponti et al. (2009) have shown in a com- 
panion paper that this emission could be well reproduced by 
a broad line from an accretion disc, while it cannot be eas- 
ily dcsc;Til)ed by scattering or emission from photo-ionized 
gas at rest. We did not investigate nor consider further 
these emission structures here, since the primary aim of 
the present work is to perform a detailed search for absorp- 
tion lines at energies above ~7 keV. We checked neverthe- 
less that different parameterizations and/or modelling of 
the complex Fe emision line structure, including the above 
weak ionized components, do not affect the results and con- 
clusions presented here. 



M. Cappi et al.: High velocity and variable iron absorption in Mrk 509 



3 



2.3. Fe Absorption Lines 

Dadina ct al. (2005) reported the detection of both rcd- 
and blue-shifted absorption lines due to highly ionized Fe 
in the BeppoSAX and XMM-Newton (non-simultaneous) 
spectra of Mrk 509. These lines were found at rest-frame 
energies of ~5.5 keV and ~8.2 keV, respectively. The first 
was seen in two BeppoSAX observations, the latter was seen 
by both satellites. Motivated by these results, we performed 
a systematic and blind search of absorption structures in 
the ~4-10 keV band for the whole XMM-Newton dataset 
available to date. 

We first fitted the 3.5-10.5 keV data with a simple 
power-law model and searched for deviations in the Ax^ , as 
shown in the first and second panels of Figure 1 for obser- 
vation 1. We then implemented a more systematic method 
to search for features in the spectrum, similar to the one 
introduced by Miniutti & Fabian (2006) and Miniutti et 
al. (2007). We fitted the 3.5-10.5 keV data with the base- 
line power-law only (panel 3) or with the power-law plus 
Gaussian emission line model (panel 4) obtained in §2.2, 
stored the resulting value, and then freezed the baseline 
model parameters to their best-fit values. We then added 
a narrow (unresolved, a — 0.01 keV) Gaussian line to the 
model, and then searched for the presence of both emission 
and absorption features by making a series of fits stepping 
the new Gaussian line energy in either the 4-10 keV or 7-10 
keV intervals and its normalization between —4 x 10~^ and 
+4 X 10^^ ph s^^ cm~^, each time storing the new x| value. 
We then derived a grid of values and made a plot of the 
contours with same Ax^ level with respect to the baseline 
model fit. Negative Ax^ values calculated in this way in- 
dicate that a better fit is reached with the inclusion of the 
line and - like in the standard steppar type of contours in 
XSPEC,- the Ax^ values can be translated in a statistical 
confidence level for the addition of 2 more parameters, as 
these do not depend on the sign but only on the absolute 
value of Ax^- We then plot the standard values of Ax^=- 
2.3, -4.61 and -9.21 relative to confidence levels of 68%, 
90% and 99%, respectively. The major difference between 
these contours and the standard ones obtained with the 
steppar command in XSPEC, is that the contours are here 
"inverted" in the sense that inner contours include higher 
significance than outer ones. These contours are shown in 
panel 3 and 4 of Figure 1, for observation 1, as well as in 
Figures 2 and 3, for all the observations. 

The probability levels shown in Figure 1 are derived 
from F-test only statistics, and do not take into account 
the number of trials used by our line blind search technique 
approach (e.g. Protassov et al. 2002). This might result in 
an overestimation of the statistical significance associated 
with each line feature. For this reason, contours derived 
with a method based on Monte Carlo simulations are shown 
for observation 1, in the inset of panel 4 of Figure 1. We 
address issues and results from this method further in §2.4 
below. 

As for the redshifted absorption structures, we noted 
the presence in the data of weak deviations between 4 and 
6 keV (see the residuals redward of the FeK emission line in 
Figure 2). After inclusion of the neutral FeK emission line 
in the fit, we found that none of these absorption features 
is, however, statistically significant with upper-limits for an 
absorption line in that energy range of EW^-15 eV (90% 
confidence). We did not confirm therefore the two positive 




5 6 7 8 9 

Energy (keV) 



Fig. 1. Panel 1 : EPIC pn spectrum of Mrk 509 (obs. 1) fitted in 
the 3.5-10.5 keV band, with data rebinned with S/N>12, fitted 
with a single power- law model; Panel 2: Ax^ residuals; Panel 
3: F-statistics confidence contours (68% (red), 90% (green) and 
99% (blue) levels) with a single power-law model with Galactic 
absorption; Panel 4'- F-statistics confidence contours (68%, 90% 
and 99% levels) after the inclusion of a Gaussian emission line 
in the model; and Panel 4, inset: Confidence contours between 
7-10 keV calculated using a Monte Carlo simulation (see text 
for details). Contours (in black) calculated with Ax^=0-5 are 
reported to indicate the continuum level. 

BeppoSAX detections of the redshifted absorption line but 
are still consistent with a scenario where the absorber is 
transient and/or patchy in nature, as proposed in Dadina 
et al. (2005). 

As for the blueshifted absorption structures, we find in- 
stead several in the data. Not only the previous XMM- 
Newton detection during observation 1 by Dadina et al. 
(2005) is confirmed, but the line is detected at a high sig- 
nificance level (see Table 3, and Figures 2 and 3) also during 
observations 4 and 6 (at E~8.2, 8.0 and 8.4 keV, respec- 
tively). In addition, there is evidence during both observa- 
tions 1 and 6 that an additional absorption line, possibly 
broader or maybe due to a blend of lines is present at even 
higher energies, at E'^9.4 keV (observed frame), i.e. ~9.7 
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Table 3. Mrk 509 absorption line parameters (rest-frame energies). 



N EW" Ax^ F-test MC" Id vjc 



obs. (kcV) (oV) (FcXXVI) 



1 


8.29 ± 0.05 
9.74 ± 0.09 


-32 ± 20 


9.6 
8 


~99% 
~98.5% 


95% 
87.5% 


Ka 
K/3 


0.172±0.005 
0.162±0.009 


2 


=8.29"= 
=9.74'= 


^ - 23 
^ - 32 












4 


8.03 ± 0.03 

= 9.74'= 


-20 ± 10 

^ - 17 


16.5 


;599.9% 


99.8% 


Ka 


0.141±0.003 


5 


=8.29'= 
=9.74'= 


^-14 
^-17 












6 


8.50 ± 0.05 


-20 ± 12 

_34+26 


10.2 
7 


~99.5% 
~98% 


95.5% 
78% 


Ka 

K/3 


0.196±0.005 

^•-^ ' -0-01 



(a) Rest-frame energy. Errors are quoted at the 90% and 68% confidence for the first and second lines, respectively. 

(b) Significance calculated through Monte Carlo simulations, (c) Values frozen at the best-fit values obtained during obs. 1. 



koV (rest-frame). Confidence contours for these absorption 
lines are shown in Figures 2 and 3. We report in Table 3 
the best-fit parameters for these lines, as well as the signif- 
icance of their detections, calculated with the F-statistics 
(column 5) and the Monte Carlo simulations (column 6). 

Remarkably, the two sets of two lines are consistent with 
the theoretical energy peaks of FeXXVI Ka (6.966 keV) and 
K/3 (8.268 keV), if blueshifted with respect to the systemic 
velocity by v-^O.lTc (obs. 1) and v~0.2c (obs. 6), respec- 
tively. The line at 8 keV during obs. 4 is consistent with 
Fe XXVI Ka at v~0.14c. If the lines were associated, in- 
stead, with FeXXV Ka (6.697 keV) and K/3 (7.880 keV), 
the requested velocity shifts would be consequently larger 
(i.e. v=0.21, 0.18 and 0.23 c for obs. 1, 4 and 6, respec- 
tively). 

We also checked the MOS results for consistency with 
the pn instrument. Remarkably, despite the considerably 
lower statistics for the MOS (due to the much lower effec- 
tive area at E>7 keV), we found evidence, albeit with lower 
significance, of absorption features in all three observations 
1, 4 and 6 and with energies and intensities consistent with 
the pn. In particular, we found statistical improvements by 
Ax^ = 13, 5 and 3, for the addition of a line at E~8.2, 8.0 
and 8.4 keV, with EW2i-67±30 eV, -24±23 eV, and -17±16 
eV (at 90% confidence for one additional d.o.f.), respec- 
tively. These results suggest these lines are not statistical 
fakes (but see also §2.4). 

The pn background spectrum is known to have two in- 
strumental emission features, one at observed energy of 
~7.5 keV due to a Ni Ka emission line and one at ~8 
keV due to a Cu Ka emission line (Katayama ct al. 2004, 
Freyberg et al. 2004). These may have a strong impact on 
the source high-energy spectrum. We checked however that 
this is not the case for Mrk 509. The source is so bright 
(F(2-10 keV)>2xl0-"erg cm^^s"!) that the background 
level was always negligibly small, even in the FeK bandpass 
and for both pn and MOS (see Figure 4), thereby excluding 
a possible artifact due to the instrumental background. 

During observations 2, 4 (for the ^9.7 keV line) and 5 
where the features were not detected, we obtained upper- 
limits on the absorption lines at energies of E~8.3 and 9.7 
keV, with EW^-15 and ^-30 eV, respectively (see Table 
3). 

Finally, there are indications that the lines do vary 
among the observations (see Figure 3, where contours are 



plotted at rest- frame energies), clearly in blueshift /energy 
(being 0.17c, 0.14c and 0.20c during obs. 1, 4 and 6 re- 
spectively) and possibly in intensity, since they were not 
detected during observations 2 and 5. In particular, albeit 
marginal, the intensity variations are found on a time-scale 
as short as ~ 100 ks, i.e. the time encompassed between obs. 
4 and 5, but one should recall the intrinsic difficulty to place 
here stringent upper-limits on the EW of lines for which the 
velocity is unknown a-priori. The line velocity variations are 
instead unambiguous (see Figure 3), on a time-scale of ~ 
6 months, i.e. the time encopassed between obs. 4 and obs. 
6. Given the mass of the black hole in Mrk 509 of Mbh ~ 
1.4 ±0.1 xlO^ M© (Peterson et al. 2004), such variabil- 
ity time-scales suggest the presence of important radial or 
transverse motions quite close to the central supermassive 
black hole as discussed below. 

2.4. On the statistical significance of the absorption lines 

The F-test detection probabilities derived here (see Table 
3) from our line blind search can be overestimates of their 
real values (e.g. Protassov et al., 2002). We therefore ap- 
plied a more rigorous method that uses extensive Monte 
Carlo simulations (see e.g. Porquet et al. 2004; Miniutti & 
Fabian 2006; Markowitz et al. 2006) in order to assess the 
detection probabilities from our line blind search. For all 
three cases in which one or two lines have been detected 
(i.e. observations 1, 4 and 6) we considered the absorbed 
power law plus narrow Gaussian emission line as the null- 
hypothesis model. 

We simulated a null-hypothesis model spectrum with 
the same exposure as the real data and subtracted the ap- 
propriate background. This spectrum was fitted again with 
the null-hypothesis model in the 3.5-10.5 keV band and the 
new parameters were recorded. Then, we generated a new 
simulated spectrum with this refined null-hypothesis model 
in order to account for the uncertainties in the model it- 
self. Therefore, we made a fit of the null-hypothesis model 
to this new simulated spectrum, considering only the 3.5 
10.5 keV band and recorded the value. We then added a 
narrow Gaussian line {a = 0.01 keV) to the model, with its 
normalization free to be positive or negative and stepped 
the centroid energy between 7 and 10 keV at intervals of 
100 eV. Each time we made a new fit and recorded the 
maximum of these simulated Ax^ ■ Then, we repeated the 
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Obsl A 






- Obs2 A 






Obs4 n 






Obs5 p 










Obs6 / 




- / 









Line Energy (keV) 

Fig. 2. 4-10 keV confidence contours for Ax^=-2.3, -4.61 and 
-9.21 (i.e., 68%, 90% and 99% levels) calculated for both pos- 
itive and negative residuals for a narrow Gaussian line in the 
data over the 4-10 keV energy range(see text for details). The 
contours at Ax^ = -1-0.5 are also reported. 



above procedure S = 10* times in order to generate a dis- 
tribution of 10'' simulated Ax^ values. If N simulated Ax^ 
values were greater or equal to the real one, then the de- 
tection confidence level corresponds to 1 — N/ S. 

The new values of line detection confidence levels based 
on the above Monte Carlo simulations are reported in Table 
3, in comparison to the F-test probabilities. We also esti- 
mated the combined Poissonian probability of detecting by 
chance the 3 absorption lines (at E'^8-8.5 keV) in 5 ob- 
servations, and this turned out to be <0.001 (i.e. a prob- 








Fig. 3. Zoomed, rest-frame 7-10 keV confidence contours for 
Ax^=-2.3, -4.61 and -9.21 (i.e., 68%, 90% and 99% levels, colors 
are as those in Fig. 1) calculated for both positive and negative 
residuals for a narrow Gaussian line in the data between 7-10 
keV (see text for details). The contours at Ax^ ~ -1-0.5 are also 
reported. 



ability >99.9% that these lines are not fake). Considering 
also the two lines at E'-~^9.7-9.8 keV, and that they turn out 
to be at the expected energies for K/3 transitions (i.e. the 
F-statistics alone should properly work for these), would 
lower the combined probability even further. For example, 
a fit with the two lines (FcXXVI Ka and K/3) with their 
energies freezed at their theoretically expected values, and 
same blueshifted velocity, yielded a Ax^ improvement of 
19 and 13, for 3 interesting parameters, corresponding to 
a probability of >99.9% and 99.5% during observations 1 
and 6, respectively. 



3. Discussion and Conclusions 

The main results of our spectral analysis are the following: 
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MOS 1—2 
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Fig. 4. pn source and background during obs. 1 (Top), rebinned 
with S/N > 12; MOS source and background during obs. 1 
(Bottom), rebinned with S/N > 6 and 3, respectively. The spec- 
tra show that the source is over an order of magnitude higher 
in counts s~^ keV~^ compared to the backgrounds in the FeK 
bandpass, in both pn and MOS. 




Energy (keV) 

Fig. 5. Best-fit photoionized models for observations n.l (top), 
4 (middle) and 6 (bottom). See text for details. 

— Wc detected several absorption features in the energy 
range ^8-10 keV in three (out of six) XMM-Newton 
observations of the briglit and himinons Seyfert 1 galaxy 
Mrk 509. The equivalent widths of the lines are in the 
range between -20 and -60 eV. 

— The most likely identification for these features being 
FeXXVI Ka and K/3, they imply (mildly) relativistic 
blueshifts in the range ~0. 14-0. 2c, and low Ka/K/3 ra- 
tios, of the order of unity. 

— There is evidence for line velocity variations, as well as 
for appearance and disappearance of the lines, on time- 
scales as short as ~100 ks. 

The above results are reminiscent of the Seyfert 1.5 
NGC1365 (Risaliti et al. 2005), in particular for i) the iden- 
tification of the absorption lines with FeXXVI Ka and K/3 
lines, ii) the small Ka/K(3 ratios, and iii) the velocity vari- 
ations detected on similarly short time-scales for a com- 
parable black hole mass. Differences are that here we only 
detect one ionization state for iron (instead of two, FeXXV 
and FeXXVI, for NGC1365), and the line velocities arc here 
about ten times larger than in NGC1365. As shown below. 



this could be due to an intrinsically more extreme outflow 
in Mrk 509 with respect to NGC1365, in terms of either ve- 
locity or ionization state, or both. Alternatively, the struc- 
ture and velocity field of radiatively-driven winds being a 
fimction of the polar angle (e.g. Proga & Kallman 2004), 
a different viewing angle as expected for their type 1 (for 
Mrk509) and type 2 (NGC1365) classification would nat- 
urally imply different average velocities along the line of 
sight. 

Outflowing material is expected to imprint a character- 
istic P-Cygni proflle in the data, where the detailed absorp- 
tion and re-emission line profile depends strongly on the 
flow/ wind geometry and density. Theoretical predictions 
based on Monte Carlo calculations have been presented in, 
e.g., Schurch & Done (2007) and Sim et al. (2008), with 
some recent attempts to fit the best available observational 
data (Sim et al. 2008, Pounds & Reeves 2008). Rather than 
attempting here such a full treament and modelling of the 
absorber with data of admittedly too low statistics for such 
a purpose, we discuss the above results in the light of sim- 
pler and more qualitative arguments. 

Following Risaliti et al. (2005) and from the curves of 
growth presented in Bianchi et al. (2005), an equivalent 
width of the Ka lines EW~20-30 eV can be obtained with 
Nh^5x 10^^ cm~^ and log^>4 erg cm s~^, with turbulent 
velocity set to zero. With larger values of turbulence, i.e. 
velocities between 100 and 1000 km/s, the column density 
needs not to be as large, with values of Nh ~5-10x10^^ 
cm~^. These values would also be consistent with the low 
Ka/K/3 ratios, which indicate that saturation effects are 
significant (see Figure 4 in Risaliti et al. 2005). Given the 
low statistics available, velocity widths as large as a few 
thousands km/s, i.e. consistent with the above estimates, 
can not be excluded here. Such high values of turbulence 
would also be reasonable for bulk motions in an outflow 
with such large velocities. 

These rough estimates were then confirmed through a 
more detailed, but still approximate, fitting of the FeK ab- 
sorption lines with two different phenomenological models: 
zxipcf. a partial covering of photo-ionised absorbing mate- 
rial that is incorporated into the XSPEC package (Reeves 
et al 2008), and XSTAR tables calculated using an input 
SED with a power-law photon index of r=2 and a turbu- 
lent velocity of 1000 km/s. We obtained good and consistent 
fits to the data. For all three observations (obs. 1, 4 and 6), 
the best-fit parameter values were log^~5±0.4 erg cm s~ , 
Nh ~2-4x1023 cm-^ and a blueshift of v-0.14-0.2c. Best- 
fit models obtained with the XSTAR photoionization model 
are shown in figure 5. In these models, the lines are pro- 
duced by only FeXXVI (Ka and K/3) and, despite the large 
columns, no other signature is expected at lower energies 
because of the extremely large ionization. 

Another possible solution, albeit slightly statistically 
worse, could also fit the EPIC data with less extreme values 
of Nh ~10^^~^^ cm~^ and log5~3 erg s~^ cm. In this case, 
the absorption lines are identified with mainly FeXXV. This 
solution does, however, require even larger velocities up to 
0.18-0. 23c. Moreover this solution appears to be slightly 
fine-tuned because a slightly lower value of ^ would im- 
ply a significant low-energy curvature and a higher value 
of $, would give rise to lines from also FeXXVI. We also 
checked but did not find any evidence in the RGS data of 
OVIII absorption lines at the above mentioned velocities, 
the data yielding quite stringent limits of EW^O.2-0.4 eV. 



M. Cappi et al.: High velocity and variable iron absorption in Mrk 509 



7 



Altogether, this solution appears to be less likely than the 
one based on FeXXVI "only" and was thus not considered 
any further here. 

The high ionization parameter, together with the varia- 
tions seen in the absorption lines can be used to infer limits 
on the distance and size of the absorber. Assuming that the 
ionization state is due to illumination from the central X- 
ray and UV source, rather than heated/illuminated locally, 
wc can use ^=iyiori/(iir^), where n is the number density. 
Lion the source continuum luminosity (integrated between 
1 and 1000 Rydberg) to estimate the maximum distance r 
of the absorber from the illuminating source. The observed 
line-of-sight absorbing column A^h is a function of the den- 
sity n of the material at a ionisation parameter ^ and the 
shell thickness Ar: Ah = nAr. Making the reasonable as- 
sumption (given the r~a dependencies of winds/outflows) 
that the thickness Ar is less than its distance from the 
source r, and combining the expression for the ionisation 
parameter, we obtain the limit r < iiora/(C Ah). Using 
Lion=2xl0''''^ erg/s, Ah=10^^ cm~^ and log^~5 erg cm s~^, 
we obtain r< 2x10^^ cm. Assuming the black hole mass of 
Mbh = 1.4 xlO^ Mq (Peterson ct al. 2004) for Mrk 509, this 
corresponds to a location for the absorber within a distance 
<500 Schwarzschild radii (r,,) from the black hole. 

Interestingly, we note that the variability of the absorp- 
tion lines is indicative of a compact absorber and rules out 
a local origin, from cither the Galaxy or the Local Group, 
for the absorber (McKcrnan ct al. 2004, 2005). In fact the 
absorber appears to be even more compact than the broad 
line region which has an estimated size of the order of ~ 80 
light-days (Kaspi et al. 2005). Analogous short term vari- 
ability of iron absorption lines velocity and intensity was 
recently found for the BAL QSO APM 08279+5255 '(Saez 
et al. 2009), suggesting that rapid changes in the geometry 
of the outflow and/or of physical parameters of the absorber 
are a fundamental characteristic of the X-ray component of 
AGN winds. 

A rough estimate of the escape velocity along the dis- 
tance of a Keplerian disk is given by the equation v'^ = 
2GM/r, which can be written as v = (rslr)^l'^c. The 
observed blueshift velocity of the absorber v~0.14-0.2c is 
therefore much larger than the escape velocity at ~ 500 
Tg, of the order of ~0.05c, as well as the escape velocity 
at 100 rg, that is ~0.1c. Assuming a constant velocity for 
the outflow and the conservation of the total mass, we can 
thus estimate the mass loss rate Mgut associated to this 
outflow as Mout ^ flnr^vnip, where v is the outflow veloc- 
ity, n is the ion number density of the absorber, mp is the 
proton mass and Q is the covering fraction of the outflow. 
Following Blustin ct al. 2005, this translates to Mout ~ 1-23 
mp^^^v, i.e. for Mrk509, to Mout ^ 2.4 x (n/iTr) wq.i 
Mq yr~^, with vq.i the flow velocity in units of 0.1 c. This 
corresponds then to a kinetic power of ^Mou^u^ ^ 6.7x 
10^'' {fl/iw) vq 1^ erg s~^ that could be injected into the 
ISM. 

Following McKcrnan et al. (2007), we can derive a sim- 
ple relation for the ratio between the mass outflow rate 
and the rate of the accretion flow, i.e. Mout /Mace ~ 6000 
(f2/47r)(uo.i/Cioo)%.ii where fl is the absorber covering 
fraction, ^loo is the ionization parameter in units of 100 
erg cm s~^ and ?7=?7o.i xO.l is the accretion efficiency of the 
black hole. For Mrk 509, this relation yields Mout/ Mace ^10 
(f2/47r), implying either a low covering fraction and/or an 



intermittent outflow to permit the growth of the black hole. 
Viceversa, even a moderately large covering factor for the 
ouflow would imply significant and energetically relevant 
mass losses into the ISM. 

The high velocity, high ionization and variability prop- 
erties of this absorber suggest that it may be associated 
with an outflow launched by an accretion disk at relatively 
small radii of less than a few hundred Scwarzshild radii. 
This adds to the long and well known evidence from the 
soft X-rays of a multi-phase warm absorber in Mrk 509 with 
logC~0.9-3.26 and Ah~0.8-5.8x lO^i cm'^ (Smith et al. 
2007,), and to the "intermediate" (logC~5.4, Ah'~5.8x1022 
cm~^ and v~0.05c) absorption component found by Ponti 
et al. (2009) based on continuum variability arguments. 
Overall, Mrk 509 shows compelling evidence for a strati- 
fied absorber, as those predicted by hydro dynamical simu- 
lations of radiativcly and/or magnetically driven outflows 
(Proga 2003, Proga 'fe Kallman 2004, King & Pounds 2003, 
Murray ct al. 1995), but with ionizations and velocities up 
to extreme and challenging values. The extreme ionization 
parameter raises the question as to how much more fully 
ionised matter is included, but undetected yet, in the out- 
flow, with its implications for the outflow mass-loss rates 
and energetics. 

Finally, the frequency in AGNs of such extreme out- 
flows is to date unclear. Proper statistical studies on larger, 
complete, samples are being performed (Tombcsi et al., in 
preparation). They are needed to secure the statistical sig- 
niflcance of these features, and to assess their frequency 
and typical characterisation in terms of density, velocities, 
ionization states and covering factor. Unfortunately current 
available data do suffer an "observational bias" against the 
detection of the highest-velocity blueshifted features in that 
orbiting X-ray telescopes are of limited spectral and sen- 
sitivity capabilities at E>7 keV. Future missions such as 
IXO will likely allow an exciting step forward in studying 
mildly rclativistic outflows in AGNs, as demonstrated by 
the first attempted simulations (Cappi et al. 2008, Tombesi 
et al. 2008). 
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